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Abstract. Results are presented for our search for warm-hot gas towards the quasar Q 0329—385. We identify ten 
O VI systems of which two are within 5000 km s“^ of Zem and a third one should be of intrinsic origin. The seven 
remaining systems have Hi column densities 10^®’^<A^(Hl)<10^® ® cm“^. At least ~l/3 of the individual Ovi 
sub-systems have temperatures T<lxl0® K and cannot originate in collisionally ionized gas. Photoionization by 
a hard UV background field reproduces well the ionic ratios for metallicities in the range 10“^ ®-10“'^’® solar, with 
possibly sub-solar N/C relative abundance. For [O/C]=0, the sizes inferred for the Ovi clouds are in some cases 
larger than the maximum extent implied by the Hubble flow. This constraint is fulhlled assuming a moderate 
overabundance of oxygen relative to carbon. For a soft UV ionizing spectrum, an overabundance of O/C is required, 
[O/C]~0.0-1.3. For a hard(soft) UV spectrum and [O/C]=0(l), the Ovi regions have overdensities p/p«10-40. 
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1. Introduction 

Numerical simulations suggest the existence of a warm-hot 
phase in the intergalactic medium, 10^<r<10^ K, which 
comprises a fraction of the baryons increasing with time. 
This phase should be mostly driven by shocks, at least at 
low redshift z (Cen & Ostriker 1999; Dave et al. 2001). 
Possible signatures of the warm-hot intergalactic medium 
(WHIM) are absorptions by high ionization species such 
as O V, O VI and O vii. These absorptions are difficult to 
detect as they either fall in the Lya forest (below the at¬ 
mospheric cut-off for z < 1.92) or in the soft X-ray range. 
Successful observations of the WHIM at low z were made 
with the FUSE, HST and Chandra satellites (e.g. Tripp 
et al. 2001; Savage et al. 2002; Nicastro et al. 2002). 

At z~2-2.5, an analysis of the Oiv/Ov ratio from 
HST stacked spectra favors a hard UV background spec¬ 
trum (thus a small break at 4 Ryd) and the inferred 
metallicity is [0/H]~—2.2 to —1.3 together with an en¬ 
hanced oxygen abundance relative to carbon (Teller et al. 
2002). Detection of individual Ovi absorbers has been re¬ 
cently reported: for systems at z~2.5 with N(Hi)~10^^ ° 
to 10^^ ° cm“^, the inferred metallicity is [0/H]~—3 to 
—2 (Carswell et al. 2002) and for N(Hi)>10^^'® cm“^ the 
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metallicity is higher, [0/H]> —1.5 (Simcoe et al. 2002). 
The main heating process of the high z WHIM is still un¬ 
clear: the more tenuous regions of the Lya forest could be 
ionized by a hard UV background spectrum, whereas the 
high column density population could be shock heated. 

A systematic, large survey of quasar absorption lines at 
high S/N and high spectral resolution is being completed 
at ESO for a sample of about 20 quasars of which half 
are at z<2.6. In this paper, we present the results of our 
search for O Vi absorbers towards one quasar of the ESO 
large programme, Q 0329—385, with several unambiguous 
cases of narrow, strong and weak Ovi absorptions. The 
observations, the selection procedure for O vi systems and 
our Ovi sample are presented in Sect. 2. The constraints 
derived from the line widths are given Sect. 3. Our mod¬ 
elling of the Ovi absorbers is presented in Sect. 4. The 
summary and conclusions are given in Sect. 5. 

2. Observations and the Ovi sample 

The quasar Q 0329—385 (zem=2.423) was observed at the 
VLT with the UVES spectrograph. The full wavelength 
coverage 3050-10400 A was obtained in two settings, using 
dichroics, with an exposure time of 6 hr per setting. The 
S/N ratio is about 30 and 100 at 3300 and 5000 A respec¬ 
tively. The resolution is 6=6.6 km s“^. A modified version 
of the ESO-UVES pipeline was used, better adapted to 
quasar spectra. A full description of the data reduction 
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Table 1. Column densities of the O vi systems in Q 0329—385 


z 

n“ 

Hi 

Ovi 

Nv 

Civ 

Si IV 

2.0615 

1 

14.26'* 

13.13 

<12.30 

<12.30 


2.0764 

1 

13.71 

13.18 

12.89 

13.23 

11.50® 

2.1470 

3 

14.70 

14.03 

<12.60 

<12.30 


2.2488 

5 

14.00'’ 

14.29 

<12.60 

12.83*’ 


2.2515 

4 / 

15.59 

14.80 

13.20 

14.43" 

12.71 

2.3139 

3 

14.23 

13.46 

<12.50 

12.45" 


2.3521 

1 

13.08 

14.04 

13.78 

13.76 


2.3638 

2 

14.84 

13.77 

12.00® 

12.53 


2.3729 

3 

15.24" 

14.25 

<12.50 

12.69" 


2.4062 

2 

14.03 

13.27 

<12.30 

<11.70 



“ Number of individual O vi components. Ions spanning 
the same velocity range but with a different number of 
components are marked: ^ —2, —1, -|-1, ® +2. 

^ Main complex. 

® 2ct detection. 


Table 2. Column densities of individual Ovi sub-systems 


z(0 vi) 

Hi 

Ovi 

Nv 

Civ 

Siiv 

|Az“| 

b 

b 

b 

b 

b 

2.24835 

12.87 

13.63 

<12.00 

12.21 


4.1x10"® 

21.8 

00 

bo 


8.9 


2.25147 

14.97 

14.37 

13.00 

14.07 

12.43 

6.0x10"® 

18.0: 

10.7 

9.5 

11.9 

16.5 

2.31395 

13.90 

13.20 

<12.30 

11.98 


0.8x10"® 

29.2 

13.4 


8.6 


2.35214*’ 

13.08 

14.04 

13.78 

13.76 


5.7x10"® 

21.8 

12.7 

9.6 

9.1 


2.36385 

14.55 

13.50 

12.00" 

12.33 


6.3x10"® 

27.5 

9.5 


10.9 



“ Difference between the redshifts of O vi and C iv. 
*’ Intrinsic system. 

2a detection. 


method will be presented in a forthcoming paper (Aracil 
et al. in preparation). The absorption lines were htted by 
multiple Voigt profiles using the VPFIT software package 
(see http://www.ast.cam.ac.uk/~rfc/vpfit.html). 

To identify the O Vi systems, we used the following cri¬ 
teria: either the two lines of the doublet are clearly seen or, 
when one line of the doublet is partly blended, the pres¬ 
ence of this line is clearly indicated by a distinct structure 
within the blend. In all cases, there are Hi Lyman lines 
associated with the O vi doublets (see Table 1). 

Our sample comprises ten Ovi systems of which two 
within Ad= 5000 km s“^ of the quasar redshift. The red- 
shift, number of components used in the fit and total col¬ 
umn densities are given in Table 1. The range of V(Hi) 
and V(Ovi) covered by these systems overlap with those 
of the samples of Carswell et al. (2002) and Simcoe et al. 
(2002). The 3a detection limits for individual sub-systems 
are typically (12, 2.0 and 1.0) xlO^^ cm“^ for Ovi, Nv 


and Civ respectively. When Civ is detected, the upper 
limits for N v are estimated using the same number of 
components and same b values as for the C iv absorp¬ 
tions. When neither N v nor C iv is detected, we assume 
6=10 km s“^ and a number of components equal to those 
of O VI. 

The individual sub-systems for which a detailed mod¬ 
elling is presented in Sect. 4 have clean O vi doublets with 
little blending which ensures a correct estimate of the b 
parameter. The velocity shift between the O vi and C iv 
doublets is always small, <5.7 km s“^, and usually larger 
than that between the Ovi and Lyman lines. This was 
already noted by Rollinde et al. (2001). The column den¬ 
sities and b values are presented in Table 2. The stacked 
velocity plots are given only for the three intervening O vi 
systems with low 6 parameters. 



Fig. 1. The complex systems at 2=2.2488 and 2.2515. The thin 
line shows the best fit model only for the lines associated with 
the two O vi main clusters; there are additional components in 
the velocity range between these two clusters. The line blended 
with Civl548 at 2=2.24835 is Mgi at 2=0.76273. 
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Fig. 2. The system at 2=2.3638 and the best ht model. 

3. Temperature 

The histogram of the 6(0 vi) values is shown in Fig. 3 for 
the 24 individual components of the ten systems presented 
in Table 1, of which five are at Ar;<5000 km s“^ (shaded 
area). Six (11) systems have line widths 6 <10(14) km s“^, 
thus a temperature T<1.0(2.0) x 10® K. In the assumption 
of pure collisional ionization, this would imply 6.3x 10“®< 
Ovi/0<2.1xlO“^ (Sutherland & Dopita 1993) and very 
high oxygen abundances. At T=(1.0, 1.5 and 2.0) x 10® 
K, the values of [0/H] are log(0vi/Hi)-|-(5.6,1.5,—0.7) 
respectively. For the systems and sub-systems given in 
Tables 1 and 2 log(0 vi/Hi) is > —1.0. Consequently, the 
six individual absorbers with T< 1.0x10® K would have 
a metallicity [0/H]>4.6, i.e. oxygen would be the more 
abundant element log(0/H)>1.5. This clearly rules out 
collisional ionization for O vi absorbers with 6<10 km s“^. 

In addition the small velocity widths constrain the size 
of the absorbers. At z~2.2, the line broadening due to the 
Hubble flow implies a maximun pathlength ^h=80 6io kpc, 
assuming Ho=65 km s“^ Mpc“^, Hm=0.3 and Ha=0.7. 
In Sect. 4.2, this value is compared to the size of the ab- 



b (km.s"’) 

Fig. 3. Histogram of the b values of individual sub-systems. 

sorbers based on the gas density derived from photoion¬ 
ization models in order to constrain the spectral energy 
distribution (SED) of the UV background spectrum. 

4. Analysis of the O vi systems 

Different SEDs were used for the UV background radia¬ 
tion at z ~ 2.2. The first one is a simple analytic fit to 
the spectrum given by Haardt & Madau (1996: hereafter 
HM96) with an energy spectral index a = 1.5 (/^ oc zz““) 
shortwards of 1 Ryd, a break of 10° ® at 4 Ryd and adopt¬ 
ing their normalization of the ionizing continuum at 1 
Ryd. We also consider the spectrum derived by Madau 
et al. (1999: hereafter MHR99) which is softer with a 
stronger break (10® "*) at 4 Ryd. Note that the analyses of 
Heii/Hi (Kriss et al. 2001), Oiv/Ov (Telfer et al. 2002) 
and Ovi/Civ (Carswell et al. 2002) suggest a hard ioniz¬ 
ing spectrum similar to that of HM96. 

4.1. The Oxygen abundance 

We used the CLOUDY code (Ferland et al. 1998) to 
derive the ionization levels of the different species and 
their column density ratios. Results are given for opti¬ 
cally thin systems, [C/H]=—1 and solar relative abun¬ 
dances. The ionization balance is very similar for any 
other values of [C/H]< —1 as the temperature, thus ionic 
ratios, is not strongly dependent on metallicity. Fig. 
4 gives V(Ovi)/V(Hi) = Ovi/Hi as a function of 
N (O vi) /N{C iv) = O vi/C iv. Variation of the metallicity 
results in a vertical shift of the theoretical curves. One sys¬ 
tem (z= 2.3521 and O vi/Hi=9.1) is clearly different from 
the rest of the absorbers: its metallicity should be about 
solar which strongly suggests an intrinsic origin. The pres¬ 
ence of strong N V associated absorption reinforces this 
conclusion. For the other systems, their positions relative 
to the theoretical curves imply metallicities in the range 
[C/H]Ri—0.5 to —2.5 for the harder ionizing spectrum. 

The column density ratios of high ions is shown in 
Fig. 5. There is a single, well defined locus for the theoret¬ 
ical models with three different SEDs (HM96, MHR99 and 
our simple analytic fit to the HM96 spectrum with a break 
of 10^ °) and solar relative abundances. Two systems, with 
detected O vi, N v and Civ doublets, lie below this locus. 
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Fig. 4. The Ovi/Hi ratio versus O vi/Civ. The red and blue 
lines show the ratios obtained with the HM96 and MHR99 
spectra respectively, [C/H] = —1.0 and solar relative abun¬ 
dances. The squares are the observed values for the ten systems 
given in Table 1 and upper limits are indicated by arrows. 



Fig. 5. The O vi/Nv ionic ratio versus that of O Vl/Civ. The 
red and blue lines refer to the same models as in Fig. 4 and, 
for the black one, the spectral slope a = 1.5 and the break at 
4 Ryd equals 10^'°. The other symbols are as in Fig. 4. 

The furthest away is the above mentioned, intrinsic sys¬ 
tem: the inferred sub-solar [0/N] abundance ratio is in¬ 
deed consistent with the elemental abundances (enhanced 
nitrogen abundance relative to solar ratios) typical of most 
intrinsic absorbers (Hamann 1997). For four out of the 
other five systems with detected Ovi and Civ doublets, 
the observed ionic ratios imply either a super-solar [0/C] 
abundance ratio, as also found by Telfer et al. (2002) for 
WHIM clouds, and/or a sub-solar [N/C] abundance ratio. 

4.2. Path length and abundances of some individual 
subsystems 

For each system the ionization parameter, U, is deter¬ 
mined from O vi/Civ adopting a HM96 spectrum and so¬ 
lar relative abundances (Model A: see Fig. 6 and Table 



Fig. 6. The Ovi/Civ (full lines) or Ovi/Nv (dotted lines) 
ionic ratio versus the ionization parameter U. Red and blue 
lines refer to the HM96 and MHR99 spectrum respectively to¬ 
gether with solar relative abundances. 


3); the pathlength is derived from the modelled phys¬ 
ical state of the gas. For the three sub-systems with 
O vi/Hi>0.1, one with A^(Hi)>lxlO^® cm“^, the inferred 
sizes are compatible with limits given by the Doppler pa¬ 
rameters (see Sect. 3) for a HM96 ionizing spectrum and 
solar relative abundances; the temperature is T=(3.0, 3.0 
and 4.8) xlO^ K for the z=2.24835, 2.25147 and 2.31395 
absorbers respectively. However, for the z=2.36385 sub¬ 
system with Ovi/Hi<0.1, the inferred size is a factor ~2 
larger than the maximun size due to Hubble expansion; 
to have //Zh< 1 implies [0/C] > 0.3. 

A second model is then investigated assuming a HM96 
spectrum and [O/C]=1.0 (Model B). The derived ioniza¬ 
tion parameters are ~3 times lower and Z<cZh for all the 
individual sub-systems. Finally, since softer spectra im¬ 
ply higher U values, thus lower densities and larger path- 
lengths, we consider a MHR99 spectrum only together 
with [O/C]=1.0 (Model C). Results for Models A and C 
are very similar as the decrease in U due to the enhanced 
[0/C] relative abundance (see Sect. 4.1) is compensated 
by the increase in U required by the larger break at 4 
Ryd in the softer spectrum. The range in metallicity of 
the WHIM absorbers presented in Table 3 is [C/H]^—0.7 
to —2.9. 

The constraints on the SED of the background ioniz¬ 
ing spectrum and the [0/C] relative abundance are thus 
linked. Even for an ionizing spectrum as hard as the 
HM96 one, we find that there must be some O vi clouds 
with a super-solar [0/C] relative abundance. Adopting the 
MHR99 soft spectrum implies [O/C]«0.0 to 1.3. 

Note that in the z=2.25147 absorber, a fraction of C iv 
could reside in the Siiv phase. Using Ovi/Nv instead to 
derive U, with [O/N]=0 and adopting the HM96 spec¬ 
trum, gives Z~l Mpc, inconsistent with the Hubble ex¬ 
pansion constraint. The latter is fulhlled for [0/N]>0.7. 
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Table 3. Abundances and pathlength of the Ovi absorbers 


2(0 vi) 

M“ 

[N/C] 

U 

H/Hi 

/(kpc)*’ 

[C/H] 

2.24835 

A 

<-0.43 

-0.50 

5.21 

6.5 

-0.67 


B 

<-0.05 

-1.05 

4.58 

0.4 

-0.84 


C 

<0.03 

-0.50 

5.21 

6.5 

-1.39 

2.25147 

A 

-0.76 

-1.08 

4.55 

47.0 

-1.13 


B 

-0.18 

-1.60 

3.98 

3.8 

-0.90 


C 

- 0.22 

-1.05 

4.58 

54.0 

-1.45 

2.31395 

A 

< 0.21 

-0.63 

5.07 

37.0 

- 2.00 


B 

<0.63 

-1.18 

4.45 

2.5 

- 2.00 


C 

< 0.68 

-0.60 

5.10 

43.0 

-2.65 

2.36385 

A 

—0.45 

- 0.66 

5.03 

141.0 

-2.31 


B 

~ 0.01 

- 1.22 

4.41 

9.3 

-2.28 


C 

~ 0.02 

-0.64 

5.05 

155.0 

-2.92 


“ Model A: a=1.5, log (Hell break)=0.5, [O/C]=0.0, 
Model B: a=1.5, log (Hen break)=0.5, [O/C]=1.0, 
Model C: Madau et al. (1999), [O/C] = 1.0 
** Hubble flow constraint: Zh<80 610 kpc. 


5. Summary and conclusions 

A sample of ten Ovi absorption systems at z=2.06-2.41 
was identified in Q 0329—385. Seven systems trace the 
warm-hot intergalactic medium, outside the close neigh¬ 
bourhood of the quasar. All the detected intervening O vi 
systems have associated H i absorption with column den¬ 
sities in the range ® cm“^. At least ~l/3 of 

the individual O vi subcomponents have low velocity dis¬ 
persions, 5<10 km s“^. These unambiguous cases of low 
temperature, r<lxl0^ K, imply a radiative ionizing pro¬ 
cess for the teneous regions of the WHIM. 

The observed Ovi/Hi and Ovi/Civ ionic ratios are 
well reproduced by photoionization models with a hard 
(HM96) spectrum and [O/C]=0, whereas the weakness 
or non detection of N v implies [O/N]>0 in most cases. 
The range in metallicity for our small sample is large, 
[C/Hjw—0.5 to —2.5, which strongly suggests a non- 
uniform enrichment of the intergalactic medium. The size 
inferred for low b absorbers of smaller Ovi/Hi ratios are 
about twice as large as that inferred from the Hubble flow 
constraint. For these absorbers, higher gas densities, thus 
sizes compatible with the Hubble flow, are obtained for 
[O/C]~0.3. 

An enhanced oxygen abundance is nearly always re¬ 
quired for a softer spectrum (i.e. a larger break at 4 ryd), 
with [O/C]«0.0-1.3, and the enrichment of the WHIM 
should then comes predominantly from massive stars. 
Adopting the MHR99 spectrum and [0/C]=I leads to a 
metallicity range [C/H]«—1.2 to —3.0. 

The A^(Hi) range of the Ovi systems in Q 0329—385 
is similar to that of the O vi sample of Carswell et al. 
(2002) and these authors also favor photoionization of the 
WHIM. They rule out a soft UV flux for a solar [0/C] 
abundance ratio, as also found for most of the O vi sys¬ 
tems in Q 0329—385 but, as discussed in Sect. 4.2, a 
softer UV flux is acceptable if [0/C] is super-solar. On 


the contrary, Simcoe et al. (2002) favor shock-heated gas. 
However, their O vi survey does not trace regions of lower 
column densities (no system with A^(H i)<l x lO^'^'^ cm“^) 
and five out of their 12 intervening O vi absorbers have 
very large V(Hi) (10^®® to 10^®'® cm“^). The Ovi sys¬ 
tems with stronger H i absorption may trace the outer 
parts of high z galatic halos where the effects of galactic 
winds may dominate the heating process (Theuns et al. 
2002), one of the models also suggested by Simcoe et al. 
( 2002 ). 

For a hard UV spectrum and [O/C]=0, the gas density 
is in the range nH=(7-23) x 10“^ cm“^, thus Ovi arises 
in regions with overdensities p/p«10-40 for flbh^Qg=0.02 
(O’Meara et al. 2001). Adopting a MHR99 soft spectrum 
and [0/C]=l leads to the same range of overdensities, 
whereas for the HM96 spectrum and [0/C]=l the over¬ 
densities are about three times larger. 

We are currently building a large O vi sample from 
the dataset of the VLT-UVES quasar large programme 
to investigate the physical state and chemical evolution 
of the Ovi absorber populations, as well as the relative 
occurence of O vi and C iv systems. 

Acknowledgements. We are grateful to R. Carswell for stimu¬ 
lating discussions and his assistance with the VPFIT s/w code. 
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